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Abstract 



o 

(N 

lf> . We show, in full generality, that the staircase method [27, 29] provides 

("^ ' integrals for mappings, and correspondences, obtained as traveling wave 

(^ , reductions of (systems of) integrable partial difference equations. We ap- 

ply the staircase method to a variety of equations, including the Korteweg- 
De Vries equation, the five-point Bruschi-Calogero-Droghei equation, the 
QD-algorithm, and the Boussinesq system. We show that, in all these 
J^ ' cases, if the staircase method provides r integrals for an n-dimensional 

?— ^ ' mapping, with 2r < n, then one can introduce q < 2r variables, which re- 

duce the dimension of the mapping from n to q. These dimension-reducing 
variables are obtained as joint invariants of fc-symmetries of the mappings. 
Our results support the idea that often the staircase method provides suf- 
ficiently many integrals for the periodic reductions of integrable lattice 
equations to be completely integrable. We also study reductions on other 
quad-graphs than the regular Z^ lattice, and we prove linear growth of the 
multi-valuedness of iterates of high-dimensional correspondences obtained 
as reductions of the QD-algorithm. 



1 Introduction 

The field of integrable partial difference equations emerged in the late nineteen 
seventies, early eighties [1, 11, 17, 23, 30]. An important, and well-studied, class 
of partial difference equations is the class of (scalar) equations that are defined 
on the elementary squares of a lattice. An example of an integrable equation in 
this class is the lattice potential Korteweg-de Vries (pKdV) equation 

(W(,m - Ui+i^,n+l)iui+l^m ~ Ui^,n+l) = a, (1) 

in which the linear terms are transformed away. Such equations are part of 
a multi-dimensional family of mutually consistent partial difference equations 
[20, 24]. A classification with respect to multi-dimensional consistency has been 
achieved recently [2, 3]. 

For lattice equations on the square, an initial value problem can be posed 
on a so called staircase: a connected path which is nondecreasing, or nonin- 
creasing. In [27] initial values are given at lattice points ui^i and mj+i.; which 
satisfy the periodicity m;.,„ — ui+p^m+p- By doing so, the partial difference 
equation (PAE) reduces to a multi-dimensional mapping. The authors of [27] 
used the linear spectral problem (Lax pair) of the lattice pKdV equation to 
derive a set of polynomial invariants for this mapping. They constructed a so 
called monodromy matrix, which is an ordered product of Lax matrices along 
the staircase over a one-period distance. This method, nowadays known as the 
staircase method, is an important tool in proving complete integrability in the 
sense of Liouville- Arnold. Here, a 2n-dimensional mapping is said to be com- 
pletely integrable if it admits n functionally independent integrals in involution 
with respect to a symplectic form [7, 42]. Thus, the number of integrals should 
be at least equal to half the dimension of the mapping. 

In [10, 21] the authors established the involutivity of the integrals for the 
mappings they introduced in [27]. Similar results have been obtained for map- 
pings derived from the lattice Gcl'fand-Dikii hierarchy [22] and for reductions 
of the time-discrete versions of the Bogoyavlensky equations [26] . In [29] more 
general staircases were given, corresponding to the s-periodicity condition:^ 

Ul^m ~ '"(+si,m+S2; (2) 

where si, — S2 G N+ arc relatively prime integers. The authors also suggested 
considering general si, — S2 G N+, see the third concluding remark in that pa- 
per. In recent work [38, 33] we have provided a unified picture for s-periodic 
reductions, with nonzero s = (si,S2) G Z x Z. In [33] it was shown how, un- 
der periodicity condition (2), any lattice equation /(m;,™, • • • ) =0 reduces to 
a system of r ordinary difference equations /(wJJ, ■ • ■) — 0, p — 0, . . . ,r — 1, 
where r is the greatest common divisor of si and S2- Also it was proved that 
the monodromy matrix, denoted C, is one of the Lax matrices for the reduction, 
that is, there exists a matrix A4 such that for (periodic) solutions of the system 



^Our notation differs from the one used in [29], where si = 22 and S2 = —21. 



the following holds, 

CnMn = MnCn+i. (3) 

In [38] a geometric description of s-reduction has been given. It was shown 
that for all s there exists a well-posed, or nearly well-posed, s-periodic initial 
value problem, for any given scalar lattice equation on some arbitrary stencil of 
lattice points. We expect something similar to hold for systems of lattice equa- 
tions, cf. [34]. Combining the two results; given the existence of a (nontrivial) 
periodic solution, after multiplying equation (3) by M.^^, we may conclude that 
the trace of C is an invariant of the mapping n i— > n + 1. In section 2 we provide 
a direct proof, in the spirit of the original work [29] , that the staircase method 
applies to any given system of lattice equations, if a Lax-pair and a (nearly) 
well-posed periodic initial value problem are known. 

Note, for equations on a square the monodromy matrix is defined on the 
staircase, and the initial conditions are given at all points of the same staircase. 
However, for equations on other stencils, and for systems of equations, the 
monodromy matrix is still given on the staircase, but the initial conditions no 
longer correspond to the points on the staircase. 

The trace of the monodromy matrix C depends on a spectral parameter, aris- 
ing from the Lax representation of the P AE. By expanding in this parameter we 
obtain a number of integrals. In relation to establishing the complete integra- 
bility of a mapping (or correspondence) obtained by periodic reduction a first 
question to ask is: does the staircase method yield sufficiently many functionally 
independent integrals? 

For the reductions we perform in section 3.1.1, of the Bruschi-Calogero- 
Droghei equation [6] 

(u;,m — Ui,m-l)(ui,m " U/-l,m-l) = (ui,m ^ U/,)n+l)(u/,)n — Uj+l.m+l), 

the number of integrals is exactly half the dimension of the mapping. For the 
one-parameter families of reductions we perform in section 3.1.2, of the QD- 

system [25] 

e;,m+l + qi+l.rn+l = qi+l,m + e/+l,m, Cl^m+iqLrn+l = qi + l,mCl,m, 

we have verified that the number of integrals is one more than half the dimension 
of the mapping. In other cases, there are fewer integrals than half the dimension 
of the mapping obtained by periodic reduction. For example, performing peri- 
odic reductions of the lattice pKdV equation (1) we find {2n -\- l)-dimensional 
mappings and (2n + 2)-dimensional mappings, for which the staircase method 
provides only n integrals, see section 4.2.1. It turns out that when the dimen- 
sion of the mapping is 2n -I- 1 it can be dimensionally reduced by 1, whereas 
when the dimension is 2n + 2 it can be reduced by 3. All n integrals survive 
the dimensional reduction and we can conclude that the dimensionally reduced 
mappings posses sufficiently many integrals for complete integrability. 

To distinguish the two kinds of reductions we say s-reduction for a periodic 
reduction of a lattice equation with period s G Z x Z to a multi-dimensional 



mapping, and we say d-reduction for a reduction of order which reduces the 
dimension of a mapping by d £ N. 

In section 4.3 we show how to pose s-periodic initial value problems for the 
Boussinesq system [22, 35] 

Wl+l^m + Vl,m — Ui,mU(+l,m, 
Wl^m+1 + Vl.m — Ui^mUl.m+1, 
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Wl^m + W/+l,m+l — Ui^mUl + l.m+l + 



Performing s-reduction with s — (n — 1, 1) wc get a 2n dimensional mapping. For 
these mappings we verified, for all n < 17, that the staircase method provides 
n — 1 integrals unless 3 divides n in which case it provides only n — 3 functionally 
independent integrals. We show, for all n, that the mapping can be 6-reduced 
if 3 divides n, and that the mapping can be 2-reduced otherwise. 

These examples suggest that if the staircase method provides r integrals for 
an n-dimensional mapping, with 2r < n, then the mapping can be d-reduced, 
with d > n — 2r. However, we do not claim the above statement is true in 
general; in examples given in [34] the staircase method gives integrals of the 
form J J' where J is a 2-integral, and it does not produce the integral J + J'. 
Recall, a function J is an fc-integral, or fc-symmetry, of a mapping if it is an 
integral, or symmetry, of the k^^ power of that mapping [13]. If one has one 
fc-integral, then one can construct k of them, or, even better, k integrals. For 
example, it is easy to see that J" — J implies that both J J' and J + J' are 
integrals. In all cases considered in this paper, the d-reduction is performed by 
introducing n — d new variables, which can be obtained as the joint invariants of 
symmetries, or fc-symmetries, of the mapping, which in turn are obtained from 
point-symmetries of the partial difference equation. This will be explained in 
section 4. 

Recently, in [4], a geometric criterion was given for the well-posedness of 
initial value problems on quad-graphs. In section 5 we will show that for 'regular' 
quad-graphs, those that permit periodic solutions, the staircase method can be 
applied. We study reductions of equation H3s=o from [2], which on a Z^-lattice 
would look like 

We will consider two different quad-graphs, namely Figure 9d and 9e in [4]. 
These quad-graphs carry more lattice parameters than the standard 7? lattice, 
and these parameters do all appear in the reduced mapping. In the second case 
the lattice parameters are interchanged by the shift on the quad-graph, and we 
find the reduction to be an alternating mapping, cf. [28]. 

For certain s-reductions the periodic solutions are given by multi-valued 
mappings, or correspondences, see [38]. The staircase method applies equally 
well in such cases, see section 4.3.1, where we perform (3,0)-reduction for the 
pKdV equation, and section 3.2.1, where we (n,0)-reducc the QD-system. Here, 



another question arises: what is the multi-valuedness of the iterates of the cor- 
respondence ? 

In general, the number of image points of the nth iterate of an ?7i-valued 
correspondence would be to". However, it has been shown that for completely 
integrable correspondences the number of images under the iterates grows poly- 
nomially, rather than exponentially [41]. In section 6 we show that for the 
correspondences obtained in sections 4.3.1 and 3.2.1 the multi-valuedness of 
their nth iterate is n + 1 and 2n, respectively. 

2 The staircase method, general theory 

Let u be a multi-component field on the square lattice Z'^ and let f [u] be a 
multi-component function of u and finitely many shifts of u. We call a lattice 
equation f [u] — integrable if it arises as the compatibility condition of two 
linear equations ip = Ltp and ip = Mip. Here, ~ denotes the horizontal shift 
I !—>■ I + I and ^ denotes the vertical shift to h^ to + 1. Thus we have 

LM = ML mod f , (4) 

which is called the Lax-equation, or zero-curvature condition. The matrices L 
and M are called Lax matrices. 

As pointed out by Calogero and Nucci [9] in the continuum case, see also 
[12], the mere existence of a Lax pair is not sufficient for integrability, the Lax 
pair has to be a good Lax pair. In the discrete setting one has to be equally 
careful, see chapter 6 in the thesis of Mike Hay [14] . From this point of view one 
could argue that the staircase method tests whether a Lax pair is 'good'. The 
Lax pair would be called good (and hence the lattice equation integrable), if it 
can be used to produce a sufficient number of integrals for periodic reductions. 

We say that a lattice equation f [u] = admits a well-posed initial value prob- 
lem if from a set of generic initial points a solution can be constructed in a unique 
way. An initial value problem is called nearly-well-posed if from a set of generic 
initial points solutions can be constructed, which can take only finitely many 
values at each lattice point. We consider periodic initial value problems. In the 
first case the solution is obtained by iterating a (finite dimensional) mapping. 
In the second case the solutions are obtained by iterating a correspondence. We 
note that if the initial value problem is well-posed the periodicity of the solution 
is implied by the periodicity of the initial values, whereas when in the case of 
nearly well-posedness the periodicity of the solutions is imposed. 

Theorem 1 Let a, s be elements of Z x Z. Suppose an integrable equation 
f [u] = allows a s-periodic initial value problem which is well-posed, or nearly- 
well-posed. Then, with C being an inversely ordered product of Lax matrices over 
a connected path, e.g. a staircase, from a to a.+ s, the trace of C^ is invariant 
under any shift on the lattice, Vi G N. 

Proof: Let a, b be two points on the lattice. Define £a,b to be the inversely 
ordered product of Lax matrices along a connected path from a to b. We have to 



show that £a,b does not depend on the path from a to b. This follows from the 
fact that every square can be passed in two ways: if L, M are the Lax matrices 
at a G Z^ (and L is a Lax matrix at a ~ a + (1,0)), then from (4) it follows 
that LM = ML = £ c is wcU-dcfincd for solutions of f ful = 0. We have 

a, a L J 

^5,b = ^b.b-^a.bi-S.a = '^b,b^a'b'^a,S 

Because the initial value problem is wcll-poscd, or nearly-well-poscd, there ex- 
ists an s-periodic solution. Now let b = a + is (i e N) , so that the value of the 
solution at a and b coincide. Then C^.h = -C*. Also, -Cj^ g = -^3,5. If we denote 

/ = Tr(£*), it is clear that wc have 1 = 1 and, similarly 1 = 1. D 

The mapping, or correspondence, which generates the s-periodic solution 
is defined by updating a set of initial values through a shift on the lattice. 
Therefore, an invariant for it is given by the trace of (an integer power of) C. 
If the Lax-matrices depend on a spectral parameter, say fc, one can expand 
the trace Tr(£') in powers of k. Each coefficient then provides an integral for 
the mapping, or for the correspondence. However, these integrals are not all 
functionally independent. 

By the Cayley-Hamilton theorem any matrix C satisfies its own characteristic 
equation P(A) = Det(A/ — £), i.e. we have P{C) = 0. Therefore, given that £ 
is a n X n matrix, it suffices to consider traces of £', with i < n. Even better, 
there are certain combinations of Tr(£*), which, generally, yield a nicer basis 
of functionally independent integrals. These are provided by the coefhcients in 
P(A) = 0. For example, if n = 2 we have 

P(A) = A2 - XTiiO + (Tr(/:)2 - Tr(/:2))/2. (5) 

Note that the coefficient of A'^ coincides with the determinant of C. 

For general n, the coefficients can be obtained using Newton's identities 

n 

ne„ = ^(-l)*"^e„_ipj, (6) 

where the power sums pk are given by pfc = .t^ + a;| + • • ■ + a;J and the elementary 
symmetric polynomials e^ are given by 

k 

I<iiil<i2 ik-l<ik ik<qj = i 

and appear as coefficients in the (Vieta) expansion 

\[{X-x,) = J2{-ire,\'^-\ (7) 

i=l 4=0 

If we denote the q eigenvalues of the matrix Chy Xi, the characteristic polyno- 
mial equals the left hand side of equation (7). Using Newton's identities (6) the 



right hand side can be expanded recursively in terms oi pk = Tr(£'^). Taking 
A; = 1, 2 in (6) we find familiar coefficients ei = pi, 62 = {Pi — p2)/2, see equation 
(5). Taking /c = 3,4 Newton's identities yield 63 = {p\ — 3piP2 + 2p3)/6, and 
64 = (pt - 6pIp2 + Spj + 8piP3 - 6pi)/24. 

For scalar equations that are defined on elementary squares, initial values are 
given on staircases. So the dimension of the initial value problem is |si| + |s2|. 
A so-called standard staircase, cf. [30, 33, 38] gives rise to a particularly simple 
mapping. In fact, any mapping, defined by a shift on the lattice, is (equivalent 
to) a certain iterate of this basic one. 

To illustrate this, we have presented 
the standard (5,3)-staircase in Figure 
1. The standard staircase is the path 
through the points between the two lines, 
the black dots. They are the points 
where initial values are given. The stan- 
dard mapping will be the shift u; „ i— > 
ui+2,m+i- Note that by this shift almost 
all black dots are shifted to another black 
dot. The one black dot which is closest 
to the dotted line, is shifted to the black 
diamond, whose value can be calculated 
using the equation on the square. 

We can also see that the mapping defined by the shift ui^m *-> ui+i,m is given as 
the third power of the standard one. To evaluate the monodromy matrix C, one 
would take the product of matrices along the same staircase on which the initial 
values are given. Since for this type of equations the matrices L, M depend 
on {u,u), (u,u), respectively, the matrix C is then automatically expressed in 
terms of the initial values. We note that one can just as well take the product 
over any other one-period long path. For example, assuming that s G N x N, 
one could consider the product 




Figure 1: (5,3)-periodic initial 
value problem for equations de- 
fined on elementary squares 



j=0 i=0 



Then, one first has to calculate the points on the corresponding path to be able 
to evaluate C. In the example given in Figure 1, we would need to calculate 
the values of the field at the white diamonds, which amounts to iterating the 
mapping 8 times. 

For equations, or systems, that are not defined on elementary squares the 
initial value problem does, in general, not lie on a staircase. Depending on the 
type of stencil and on the particular periodicity condition, there could either be 
more, or less than |si| -|- |s2| initial values. In [38] it was shown how to write 
down, for a given scalar equation on an arbitrary stencil, a piece-wise linear 
expression (as a function of s) for the dimension of an s-pcriodic initial value 
problem. 



Note that the monodromy matrix is still a product of |si| + |s2| matrices (if 
the product is taken over a staircase, which is the sensible thing to do). Also 
note that, in general, the Lax matrices depend on u and a number of shifts 
of u. Therefore, for certain choices of s, one needs to determine a number of 
points, by iterations of the mapping, or correspondence, in order to evaluate the 
monodromy matrix in terms of the initial values. It might be possible to avoid 
this by using the equation to change the [u]-depcndence of the Lax-matrices. 
However, that would have to be adjusted to the particular choice of s. 

3 The staircase method, appUcations 

In this section we will apply the staircase method to a variety of different equa- 
tions and different reductions. Firstly, we consider a 5-point equation, i.e. an 
equation which is not defined on the square. Second, we calculate integrals for 
reductions of a 2-component system of equations. We first perform reductions 
which give rise to mappings and then we also perform reductions that yield 
correspondences. 

3.1 Mappings with a sufficient number of integrals 
3.1.1 The Bruschi-Calogero-Droghei equation 

In [6] one finds the five-point equation E{u, u, u, u, u) ~ 0, 

where ^ 

E — (u — u){u — u) — (u — u){u — u) , (8) 

Figure 2 see [6, Equation (4a)] in which we have set a'-'^^ — 0. The 

equation is defined on the stencil given in Figure 2. 
The two recursive formulas [6, Equations (la), (6)], with coefficients [6, Equa- 
tions (5), (7)], yield the following Lax-pair 




L{u,u,u,u) = 



k + u~u {u — u){u — u) 



1 



/ 1 (u-u) \ 

M{u,u,u) - ( (m-m)-i l-{k + u-u){u-u)-^ y 

Here we have denoted the spectral parameter {x in [6]) by fc, which we do 
throughout this paper. The method laid out in [38] tells us how to pose well- 
defined s-periodic initial value problems for this equation. This can be done 
for all s = (si,S2) such that 52(^2 — 2si) 7^ 0. The dimension of the periodic 
solutions is given by the following piecewise- linear function 2max{|s2 ^sij, |si|}. 
We apply the staircase method to a few reductions, in the different regions 
distinguished by this function, see [38, Figure 10]. 



(0,3)-reduction 




Xi X2 



We assigning initial values as in Figure 3. They are in- 
dicated by the black dots. We update them using the 
right-shift, the values of 52,a?4,X6 can be determined us- 
ing the equations indicated by the zig-zags. We get a six- 
dimensional mapping 

X2Xe+X3X4—X3Xe—X4Xe 
x2—xe 

X2X5+X2Xe—X2X4,—X5Xe 




X2—X4 

XiX4-\-X2X4^ — XiX2—X4XG 



Figure 3 



The monodromy matrix, which we take from X2 upwards 

to X2, is £ = M {xq, X2, X5)M {X4, Xq, X3)M {X2, X4, Xi) . 

Three functionally independent integrals for this mapping can be obtained from 
the coefficients in its characteristic polynomial (5). They are (xi — 0:4) (2:2 — 
2:5) (2:3-2:6), {xi-xe){x2-X3){x4,-X'i), and (xi-a;5)(2:2-a;4)-|-(a;3-a;5)(a:4-a:6). 



(l,3)-reduction 



Xl X4 



Assigning initial values as in Figure 4, and updating them 
using the down-shift, we get a four-dimensional mapping 




(X1,X2,X3,X4) 1~^ (X2,X3,X4 



X1X4 -1- X2X3 — X1X3 - X3X4 , 
X2 - X3 



Figure 4 



(9) 

In this case we have to first calculate a few values of 

the field at points close to the staircase in order to 

evaluate the monodromy matrix. We calculate Xi by 

solving i<^(x2, Xl, X3, X4, Xl) = 0, and we find xi from 

£;(xi,xr,X2,X3,xT) = 0. 

Two functionally independent integrals are obtained from the coefficients in 

(5) with C = M(x2,xi,xi)M(x3,X2,xi)M(x4,X3,xi)L(xi,X4,X2,xi). They are 

Ji = (xi - X2)(x3 - X4), and J2 = (xi - X3)(x2 - X4)/(x2 - X3). 



(2,3)-reduction 




Figure 5 



We assign initial values as in Figure 5 and update 
them using the diagonal shift u M- w. The val- 
ues X2, which equals X4, and X3 are determined by 

£^(X3,Xi,X2,X4,X2) = 0, and i?(x4,X2,X3,5f2,X3) = 0, 

successively. We find the same four-dimensional map- 
ping as in the previous case. The monodromy matrix 

M(X4, X2, Xi)M(x3, X4, X3)L(X3, X3, X2, X2)M(x2, X3, X2)- 

L{x2,X2,X4,xi) yields the same integrals Ji, J2. 



(2,-l)-reduction 




Assigning initial values as in Figure 6, and updating them 
using the right-shift, we get a six-dimensional mapping 

Xi M^ Xi+i, i e {1,2,. . .5}, 

xiXi + a;22;4 — a;ia;2 — x^xq 



xi X2 xa 



Xl 



Xq ^ 



X^ — Xq 



Figure 6 



The trace of the monodromy matrix 

M"i(x4,X6,X3)i(x5,X6,X3,X2)L(x4,X5,X2,xi) yields 

two functionally independent integrals 
These are (xi — X6)(x2 — X4)(x3 — X5) and (x2 — X6)(x3 — X5) + (xi — X4)(x2 — 
X4)-|-(x3 — X4)(x4 — X5. A third functionally independent integral, (xi — X4)(x2 — 
X4)(x2 — X5)(x3 — X5)(x3 — Xg), is obtained by taking the determinant of the 
monodromy matrix. In the previous cases, the determinant does not provide a 
functionally independent integral. 



3.1.2 The QD-algorithm 

The quotient-difference (QD) algorithm. 




e + q = q + e, 
eq = eq, 



(11a) 
(lib) 



Figure 7: 
QD-type system 



is used to construct continued fractions whose convergents 
form ordered sequences in a normal Fade table [8], and to 
find the zeros of a polynomial [15]. 
It is also called the time-discrete Toda molecule [19]. It is an integrable two- 
component equation defined on the stencils depicted in Figure 7, where we 
associate two values, e on the left and and q on the right, to every point on the 
lattice. 

A Lax-pair for the QD-algorithm can be obtained from relations between so 
called higher adjacent orthogonal polynomials [5], cf. [25, equations (3,4)] 



kP = P + qP, P = P + eP. 
With ** = (P, P) wc have $ = L* and $ = fc^^M*, where 



L{e,q) = 



-q k 
-q k — e 



M{e,q,q) = 



k 

—q k — e + q 



This (small) Lax-pair differs from the (big) Lax-pair obtained in [25, equation 
(9)]. A big Lax-pair incorporates a particular choice of initial values and we 
would like to consider general periodic initial value problems. 
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In [38], it has been shown that for ah s = 
(si, S2) G Z X Z such that 52(52 — 2si) 7^ 0, there 
exists a weh-posed s-periodic initial value prob- 
lem, with dimension 

2max(|si + S2I, |si|). 

This function tells us there are two different re- 
gions, where the dimension is given by a different 
linear function of the periods (up to a sign) . How- 
ever, to pose the initial value problems one has to 
distinguish four different regions, as depicted in 
Figure 8. 
We note that for all s initial values can be given on (part of) a standard staircase. 
In the Figures that follow this will be indicated by a dotted line. We will present 
four examples of families of periodic reductions, where the dimension depends 
on an arbitrary variable n G N, one family in each of the different regions. 




Figure 8: Distinct regions 
for the QD-stencil 



(0, n)-reduction 



X2n-1 



a;2n-3 




Figure 9: Vertical period- 
icity. 



We take initial values (e, q)o,i = {x2i-i, X2i)^ with 
Xk — Xm ii k = m mod 2n, see Figure 9. 
Updating these values to the right gives a 2n- 
dimensional volume-preserving mapping, with 
i = l,2,...,2n. 



3^21-^1 ^ X2i+1 



X2i+3X2i+i X2i+lX2i+2 



X2i+1 



X2i-1 



X2i 



X2i+lX2i+2 
X2i-\ 



Here we first used equation (llfo) to find the im- 
ages X2i^ and then equation (11a) to find X2i^\. 
The nionodromy matrix is 

C ^M{x2n~l,X2n,X2)M{x2„^3,X2n~2,X2n) ■■■ 

■ ■ ■ M{x3,X4,xe)M{xi,X2,x^). 



Wc have verified up to n = 9 that the coefficients in the fc-expansions of its trace 
and determinant yield n+1 functionally independent integrals. The lowest non- 
trivial mapping, taking n = 2, is 



{xi,X2,X3,X4) H^ .T3 + 



a;ia;2 2:3X4 X3X4 



X3 



xi 



x\ 



,Xl 



0:3X4 X\X2 X\X2 



X\ 



XZ 



xz 



which admits the three functionally independent integrals x\ + X3, X2X4, x\X2 + 

X3X4 - XxX-i- 
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(3n, 2n)-reduction 

We pose initial values as in Figure 10 and update then using a horse-jump 
{I, to) i-^(/ + 2,to — 1). Note that one should first update X6i+2 before a;6i+4- 

Thus we get a 6n-diniensional map- 



•^6n — 3 ^6n — 2 



6n- 


5a;6n- 


-4 i 












oo 


OD 


•0- - ■ 


y^n Z 


CO 


CO 


CO 


CO 


cx> 


Xl X2 


.J^i 


«> 


CO 


OD 


CO 


CO CO 


X5 


- -G* 



Xq 



ping, with i ~ 0,1, . . . ,n — 1, 

2^61+72^61+9 

a;6j+2 1-^ , 

2^61+6 
2^61+3 '^ 2;6i+6, 



X7 



Figure 10: Initial values on part of 
a standard (3n, — 2n)-staircase, and the 
image points y, z of 2:611-4, 2;6n-2- 



2;6i+4 1^ 2:61+8 + 2:5^+9(1 - 

2^64+5 '^ ^6i+Sj 
Xei+6 '^ 2:6i+io, 



X6i+7 , 

XQi+e ' 



where the subscript on x is taken 

modulo 6n. 
This mapping is measure-preserving with density (nr=i ^61-3)^^- We obtain 
integrals by expanding the trace and determinant of the nionodroniy matrix 



Jl M ^(2:61+7, 2:6j+5 +.2:6i+6 - •2:6i+7,2:6j+6)-^(2:6i+5, 2:54+2 +2:54+4 -2:61+5) 
■M" (2;6i+5 , 2;6i+2 + 2:6i+4 — 2:6i+5,a:;6j+4)i(2:6i+2,2:6i+3) 

■LyX^i^X, XQi^i + Xgi — Xgi^i) 

in powers of the spectral parameter k. We verified that, up to n = 3, 3n + 1 of 
them are functionally independent. For n = 1 the mapping reads 

/ \ f ^12^3 , /I 2:1 \ 

(0:1,0:2, 2:3, 2:4, 2:5,2:6) H+ 2:5, ,X6,X2 +2:3(1 ), 2:2, 2:4 , 

V 2:6 2:6 / 

which admits the following four functionally independent integrals 

2:6 + 2:4 + 2:2 + 2:3,^3 (2:4 -2:5) (2:1 -X6),X3(x6+X4-Xi)+X6(x4 + 2:2 -2:5), X1X2X3X5. 
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{2n, — 3n)-reduction 



90- -■ 

Xin- 


XAn 
-3 : 


-2 

00 CO 


OO 


ccIh-I 


- -GO CO 

1 : 


00 


CO 


\ 


F4"co 


CO 


OD 


• 
X 


o...\-^2 

1 : 


00 


CO 





*> 

Xz 



Wc choose initial values as in Figure 11, 

e2i+l,-3i-l = a;4«+3j <l2i+2,-Zi-2 = X4i+4, 

with Xk — Xm a k = m mod An. They are updated 
by shifting {l,m) i— > (Z + l,m — 1). This yields the 
4n-diniensional mapping 



CO CO OD OD 



■X4 



X4i+1 >-> X4i+3, 

Xii+^Xii+Q 



X5 

Figure 11: Part of the 
(2n, — 3n)-staircase 

The monodromy matrix is 



X4i+2 



2^44+3 



X4i^4 
Xii+bXii+e 



+ Xii+i — a;4i+3, 



Xii+i 
X4i+4 ^^ X4i+Q, 



which is measure-preserving with density Y[i=i ^'H- 



Y\^M ^{x4i+5,Pi,X4i+4)M '^{ei,X4i+4,qi)L{x4i+s,ri) 

4=0 

where e^ = X4i+5X4i+6/X4i+4, Qi = e^ +^4^+4 -X4j+3, n = X4i+i + X4i+2 - X4i+3, 

Pi = Zi + X4i+4 — X4i+5, and Zi — X4i^^ri/x4i+4. We verified up to n = 3 that 
its trace and determinant yield 2n + 1 functionally independent integrals. For 
n = \ the measure preserving 4-dimensional mapping reads 



XiXZ X1X2 

(xi,X2,X3,X4) h-> 2:3, |-a:;4 - x-^, ,X2 

' 2:4 X4 



which admits the following three functionally independent integrals, 

xix^ {xi - X4){x2 - x^){x:i - X4) 



2^2 +a;4 - 2:3, 



(11) 



X4 



Xi 
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(1, 1 n)-reduction 



2:3 

Xi 



Xn+2 
I 00 

Xn+1 




lX2n 
CO Xn^2n—1 
Xn-l'-X2n-2 



We choose initial values as in Figure 12, with j = 1,2, . . . ,n, 
and i € 7j, 

Updating by the up-shift yields a 2n-dinicnsional mapping 
Xii~>Xi+i, i e {1,2,. . .,2n - l},i 7^ n 

, X2Xn+l 
X„ ^ Xi H X2„, 



Xi 



a;2n >-> 



2^22; Tl+l 
Xi 



Figure 12 
The monodromy matrix is 



which is measure-preserving with density xi. 



M (Xi,(x„ +X2„_1 - X2„) ,X„+i)(TTM (Xi,X„+i_l,X„+i)) 

Xi -'■-'- 

^ i=2 

.,_i. , a;2X„+i X2X„+i X2„ 

• M (Xi H X2„, X2„, 2 )^[Xl, [Xn + X2n-1 - X2„) ), 



Xi 



Xi 



Xi 



whose trace and determinant yield n + 1 functionally independent integrals, 
which we verified up to n = 8. For n ~ 2 we find mapping (11) again, under 
the change of variables 

(xi,X2,X3,X4) 1~> (X4,X2,X1,X3). 

3.2 Correspondences with a sufficient number of integrals 

There are certain lines in the (si, S2)-plane where a periodic reduction yields a 
correspondence instead of a mapping, see [38] . Here we impose the solution to 
be periodic, which is not implied by the periodicity of the initial conditions. For 
the QD-system, we find correspondences on lines given by 52(52 + 2si) — 0, cf. 
[38]. 

3.2.1 The QD-algorithm, (n, 0)-reduction 

We consider horizontal staircases for the QD-system. Solving the non-local (or, 
implicit) scheme we find rational expressions for two-valued correspondences. 
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_C2^2n ^^^2 


CO 


C2n-1 Xi 


X3 


X2n X2 


-..•.^-4 


a;2n-i Xi 


Xi 



As initial conditions wc take 



^rn,0 -^2771— 1 7 Qm,0 X2'tm 



where the index on x is taken modulo 2n, 

see Figure 13. They are updated by the 

Figure 13: Initial values and up-shift. We assume the image is periodic 

their images under the up-shift. ^^^^ ^^"^ ^^me period as the initial values, 



i.e. we also take the index on x modulo 2 



n. 



We have the following 2n equations for the 2n unknowns Xi (the reader might 
like to draw a few of them into Figure 13) 

x^ix^i = a;2i-ia;2i+2, x^i + x^2 = X2i+i + a;2j+2- (12) 

We first solve for the odd variables, thereby obtaining a set of n equations for 
the even variables 

^ X2i-lX2i+2 , ^ ^ 

a;2j-i = ^^ = a;2i+i + a;2i+2 - a;2i+2- 

a;2j 

We write f^^ = m2i+2ix2i), where 

mk{z) = Xk-i + Xk . 

z 

Now 'x2k must be one of the two fixed points of the Mobius transformation 

Mk = m-2fcm2(fc-l)TO2(fe_2) • • •TO2(*;-n+l) 

= m2fe • • • m2(fe+3)m2(A;+2)m2(fe+i). 

The first fixed point of Mk is given by a;2fc-i, as vai{xi-i) = Xi-i. This gives 
us one way of updating our initial values, i.e. the linear map a2n '■ 

X2^-\ y^ X2^+2, , , 

X2i l-> X2t-\, 

taking i — 1, 2, . . . , n, assuming the index on x to be periodic modulo 2n. 
The other fixed point of M2k is Zfc '■— X2(k+i)Q2k/Q2k-2^ where 

n-l A-l \ /n-l \ 

Qk{x) = Y. [li ^2.+2+/c n ^2^ + 1+'^ • (14) 

This follows from 
Lemma 2 



m2kizk-i) = Zk 
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Proof: In terms of Q the statement is 

X2k-zQ2k-A + X2k+2Q2k = {x2k-l + X2k)Q2k-2 (15) 

From the definition (14) it follows that 

n n 

X2n-lQ'^'^ + \{X2J - Qo = a;4gr' + n^2j-l. (16) 

Therefore we have 

X2n-lQ'-2 +XaQ2 

n—1 n+1 

= X2„-i(x2(3o"^ + n ^2j-i) + xa{xiQ'^-'^ + ]J a;2j) 
== X2(x2„-iQo"^ + W X2j) + xiixiQ'^^'^ + Y[ a;2j-i) 

3=2 3=2 

= (Xi +X2)Qq, 

which is equation (15) with fc = 1. This implies that equation (15) holds for all 
k, as we may shift Xi t-^ x^^2{k-i)- D 

The mapping that corresponds to the fixed point Zk will be denoted T2n '■ 

Q2i-2 

(17) 



X2i-1 H> X2i-1- 

'^2^ 



^2i 
X2i l-> X2i+2-7Z;i • 

'^2i-2 

Since a Mobius transformation has at most two fixed points we obtained a two- 
valued correspondence (cr, r). Integrals for this correspondence are given by the 
coefficients of the fc-expansions of the trace and determinant of the monodromy 
matrix 

L{x2n,X2n-l) ' • ' L{x4:, X3)L{x2,Xi). 

For all n < 10 we found n+1 functionally independent integrals. The lowest 
non-trivial case is n = 2. Explicitly, both mappings a^ : 

(xi,X2,X3,X4) H^ {X4,XI,X2,X3) 

and T4 : 

, , f X3 + X4 a:;i + a;2 xi -|- 2:2 2:3-1- X4 

(Xl,X2,X3,X4) ^ [Xl ; ,X4 ; ,X3 ■ ,X2 ■ 

\ Xl+ X2 X^ + X4 X3+ X4 Xl + X2 , 

admit the three invariants xi -|- X2 -|- X3 + X4, X1X3 + X2X4, and X1X2X3X4. 
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4 Reduction of order 

In all cases we have encountered so far, the staircase method provided a sufficient 
number of integrals, that is, at least half the dimension of the mapping. Still, 
the dimension of those mappings may be reduced. This can be done using a 
symmetry of the lattice equation. In certain cases the number of invariants 
remains the same, whereas in other cases it drops. 

For reductions of other equations, which we encounter in this section, the 
number of integrals provided by the staircase method is not sufficient. However, 
in these cases, after reduction of order the number of integrals will suffice. We 
observe that the number of dimensions to be reduced varies with the period s. 
This can be understood exploiting symmetries of the PAE that give rise, for 
certain periods, to /c-symmetries of the mappings. 

4.1 Mappings with sufficiently many integrals, revisited 

4.1.1 The Bruschi-Calogero-Droghei equation 

For the mappings obtained in section 3.1.1 the number of functionally indepen- 
dent integrals is exactly half the dimension of the mapping, sufficiently many 
for complete integrability. We note that the equation E = 0, cf. (8), admits 
two Lie-point symmetries u i-^ u + e and u t-^ Xu. This yields two symmetries 
for the mappings, which can be used to reduce the dimension of the mappings 
by 2. The integrals we have given only admit the first (translation) symme- 
try. However, certain homogeneous combinations of them also admit the second 
(scaling) symmetry. Therefore, applying 2-reduction to the examples in section 
3.1.1 produces 2-, respectively 4-dimensional mappings with 1, respectively 2 
integrals. For instance, using reduced variables 

Xi - X2 X2 - X3 

Zl = , 22 = , 

X2 - X3 X3- X4 

the 4-dimensional mapping (9) reduces to 

{Z1,Z2) ^ (z2, — ], (18) 



which has one integral 



Ji (Z1 + 1)2(Z2 + 1)2 



Jl Z1Z2 

Note that the 4th iterate of (18) equals the identity and hence the reduction 
provides an explicit solution for mapping (9). If the (n — l)st iterate of the 
mapping (9) is denoted (x„,x„+i,x„+2,a;„+3), then 

Xn = Xl+ 1——\{X2 - Xi) + [——\{X3 -X2) + [j\{x4- Xi) 

. n-1 {x2 -a:i)(a;4 - X3) 

L ^ J ) 

4 xs - X2 
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where [ J denotes the floor function, sending x to the largest integer below x. 
This solution can be obtained similarly to the solution of the (3,l)-reduction of 
lattice pKdV given in the appendix of [39] . 

4.1.2 The QD-algorithm 

As the QD- system (11) admits 1 (scaling) symmetry, all mappings can be 1- 
reduced. The reduced mapping has one integral less, as only homogeneous 
combinations with scaling eigenvalue are invariant under scaling. For example, 
the mapping (11) with new variables yi — Xi/x4, i — 1,2,3, reduces to the 3- 
dimensional mapping 



/ ^ , / 2/3 1-2/3 
(2/1,^2,2/3) f-^ — ,2/1 H ,2/1 , 

\2/2 2/2 



(19) 



which has invariants 



h = 



2/12/3 



2/2-2/3 + 1 



l2 = 



(2/2 -2/3) (2/1 - l)(2/3 - 1) 

(2y2 - 2/3 + 1)2 



This mapping has a lot of periodic points. The first few are given in Table 
1, where a,b are free parameters. Also, one can show that the orbit of {a,h,b) 
converges to the periodic orbit of (l,a6, 1), which has length 3. 



orbit length periodic points 



1 
2 
3 

4 
5 



(0,-1,0), (a, 1, a) 
(0,1,2), (2,-1,0), (0, a, 0) 
(l,a,l), (l,a,a), (a, 1,1) 
(a, 6 -1,6) 
(a,5(6+l)/(a + 6-l),6) 



Table 1: Periodic points with orbit length smaller than six. 

Introducing variables a; = yi , y = j/3 on a level set of the first invariant Ii = z 
yields the mapping 



T ■■{x,y)^ { 



zy 



xy + zy — z 



This map can be written as a composition of two involutions, namely t = zooii, 
where 

lo : (x,2/) h^ (y,x), h: {x,y) ^ { ^ -^,y 



h - hx 



with / = Av X Bv, V = {y'^,y, 1)*, which has invariant 



w ■ Av I2 
w ■ Bv z ' 
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-1 


0\ 




(I 








-1 


1-z 


z , 


B = 














z 


—zl 




\o 









where w = (a;^, x, 1). Thus, it is a special case of the 18-parameter QRT-family 
of planar maps [31, 32], with 



A 



Reduction of order also works for correspondences. In reduced coordinates 
iji — Xi^i/xi, i ~ 1,2, 3, the four dimensional mappings (T4, T4 reduce to 

<j: {yi,y2,y3)-^(-A,-) (20) 

\y3 ya ysj 



and 



t: (yi,2/2,y3 W J/37 ] ^'2/27 — \ — r^,yi 

V {y3 + y2r (^3 + 2/2)^ 



(21) 



which admit the invariants 

2/2 + yiys yiy2yz 



(1 + yi + 2/2 + 2/3)2 ' (1 + y^ + y^ + y^Y ' 

4.2 Mappings with insufficiently many integrals 

We will next encounter reductions whose dimension is greater than twice the 
number of functionally independent integrals provided by the staircase method. 

4.2.1 The potential Korteweg-de Vries equation 

The matrices 

u —k — uu \ f u —a — fc — uu 
1 -u ) ' \l -u 

form a Lax pair for the lattice pKdV equation 

{u-u){u-u) = a. (22) 

4.2.2 (n — 1, l)-reduction 

xi X2 x^ We consider the initial value problem 



Xn—l 



Xn 



Xn '"(ri-l)fc+i-l,fc 



Figure 14 with i = 1, . . . ,n and fc £ Z, see Figure 14. We 

find, for all fc, U(„-i)fc+„,fc ^ xi + a/{xn - X2). 
Hence, the right-shift induces an n-dimensional mapping, 

(xi,X2,...,X„) ^ X2,X3,...,X„,Xi H I , (23) 

V Xn- X2 I 
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which is volume-preserving when n is odd and anti- volume-preserving when n 
is even. 

We have verified up to n = 17 that the trace of the monodromy matrix 
provides \_{n — 1)/2J functionally independent integrals. So, for the odd dimen- 
sional mappings we need one reduction, but for the even dimensional mappings 
we need two. In fact, as we shall see, in the even case there exist three. How to 
explain this? As we will see below, the KdV-equation has a Lie-point symmetry 
that does not depend on the lattice variables. This symmetry gives rise to one 
reduction for both the odd and the even dimensional mappings. Also, there 
are two symmetries that do depend on the lattice variables. These yield two 
2-symmetries of the mapping if and only if its dimension is even, giving us two 
more reductions. 

Equation (22) has the following symmetry u i— >■ u+e. This yields a symmetry 
for the mapping (23), whose infinitesimal generator is 

z— 1 

Let jji = Xi — Xi+i for i = l,2, ...,n — 1. The functions y, are annihilated by the 
vector field (24) , they form a set of n — 1 functionally independent invariants of 
the symmetry generated by this vector field. Taking the y as new variables the 
mapping (23) reduces to 

(yi,2/2,---,yn-i) ^ Iy2,y3,---,yn-i,-^yi+ _i • (25) 

In addition, equation (1) also has the following two symmetries 



Suppose now that n is even. Then the above symmetries of the PAE (22) give 
rise to 2-symmetries of the mapping (25), with generators 

It can be verified that of n — 3 functionally independent joint invariants of the 
above three vector fields (24, 26) are 

Qi = {xi - Xi+2){xi+i - Xi+s), I = 1, . . . , n - 3. 

We will take the qi as reduced variables and perform the reduction. Let us 
define, with fc, to S N, 

Tpm ^1 + fe ~ X2rn+3+k 
^k '■- • 



X2m+l+k — X2m+3+k 
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The F™ satisfy the recurrence relation 

(F™ - l)q2m+k = '?2m+fe-l-F" , (27) 

with initial condition Fj^ = 1. Therefore the FJT can be expressed in the qi, 
with i < n — 3, when m < n/2 — 1 and k < n/2 — 2m + 3. We have the following 
equalities 

Xi = a;„_i H 1'-_^ , J odd, 

-^n — 2 '^n 

Xi = x„_2 H (a;„_2 - a;n)-f'i_i , i even, 

which is an inverse reduction. In terms of reduced variables the mapping is 
Qi^qi+i, i e {l,2,...,n-4}, 

qn-3 ^^ —qn-3-fo +a/r-^ . (ZS) 

An explicit expression for F™ in terms of qi is^ 

m+l 2m+l 

^r = E n '^i-:!^.' (29) 

i=l j=2i 

since this expression solves the recurrence (27) with Fj^ = 1. The mapping (28) 

is anti-measure preserving with density Y[i=i 92i- 

At n = 4 the reduced mapping is qi h-> a — qi, which admits one integral, 
qi{a — qi). The second iterate of this mapping equals the identity. Note that 
this enables one to exphcitly solve equation (23) with n = 4, cf. [39]. What 
happened? Well, the joint invariant qi turns out to be a 2-integral of the 
mapping. Let us define another set of functions 

2m+l 
1=1 

They can be expressed in terms of the gi, with « < n — 3, if to < n/2 and 
k < n/2 + 4 — 2m, using the recurrence 

Hr^H^r'+q2m-i+kFr\ (30) 

with initial condition H^ — 0. For all n the n-dimensional mapping J„, see (23), 
admits the 2-integral Hq . If the image of x under (23) is denoted x, then 

rii^^/2 _ ^ _ TT(n-2)/2 



^The function _F™ is closely related to the function Q given by equation (14), we have 



^r - 1 = Qr+"i'(9)(n ■=! 92,+™)-^ 
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An explicit expression for H^ in terms of the qi is 

m m+1 — « 2«— 2 



Hr-T.T. n4Xi..> 



4=1 j=l 1 = 

since this solves the initial value problem (30). At n = 6 the 3-reduced mapping 
is 

[qi,q2,q3) ^ (q2,q3, \-a — ■ — ). (31) 

92 92 + 93 

The staircase method provides two integrals 

f , , 9lfe + 93)^ . , , 9193 w , , 9193 X 
9193(92 + 93 - a H ), (91 + 93 H )(9i + 93 - " H ), 

92 92 92 

of which the latter can be expressed in terms of the 2-integral Hq = 91 + 
93 + 9193/92- We can take the 2-integral as a variable. In terms of 91,92 and 
p — Hq mapping (31) becomes 

(9i,92,p) ^ (92,92(p-9i)/(9i +92),a-p) 
which has integrals 

\ (p-9i)(p + 92 -a),p{a-p). 

91 +92 

In general, with n > 2 even, the 3-reduccd mapping can be written as 

p _ ^("-4)/2 
(9l,92,...,9n-4,p) 1-^ (92,93,---,9n-4, (TT^^yy^— ' " ^ -P) ' 

where we used (30) to solve p ~ Hq for 9„_3. 

We note that all [(n — 1)/2J functionally independent integrals we have 
calculated for n < 17 dimensional mappings survive these reductions. Taking 
n = 2m + 1 odd, the reduced mapping is 2m dimensional and has m integrals. 
With n — 2m + 2, the reduced mapping is 2m — 1 dimensional and has m 
integrals. 
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4.3 The Boussinesq system 



u V u w 
» «> 



u w 

«> 



U U V 

w m 



-* 



U V 



^ 




(a) 



ib) 



(c) 



(d) 



The Boussinesq system [22, 35] 



w ~ uu — V 
w = uu — V 



w = uu 



V - 



u — u 



(34a) 
(34b) 

(34c) 



is defined on the square as depicted in Figure 15(a,b,c). In 
Figure 15(d) we have depicted the consequence of equations 
(34a,34b), 

uu ^v — uu — V. (34d) 

From the Boussinesq system one can eliminate the variables 
V and w using the identity 

iv + w) — {y + w) — {v — v) + {w — w) 



Figure 15 



to get a 9-point scalar equation on a 2 x 2 square, called 
the Boussinesq equation, cf. [22, equation 1.3]. 
We denote u — {u,v,w). A Lax-pair for the system (34) is given by, cf. [35], 



L, 



M„ 





( 


—u 


i 


U 


\ 


: L{u, w, u, v) = 




~v 





1 


, 




y uw 


-vu — k 


— w 


u 


/ 




( 


—u 




1 





M{u,w,u,v) — 




— V 







1 




\ uw 


-vu + j- 


-k 


—w 


u 



(33) 



Invariants for traveling wave reductions of the system can be obtained by ex- 
panding traces of powers of the monodromy matrix. Since £ is a 3 x 3 matrix a 
full set of functionally independent integrals can be obtained from fc-expansions 
of the coefficients in 



Det{XI - C) =A3 - X^TriC) + A 



Tr(/:)2 - Tr(/:2 



TrjCf - 3Tt{C)Tt{C^) + 2Tt{C^ 
6 



(34) 



cf. section 2, in particular equation (7). However, due to the fact that both 
Lax-matrices have a constant determinant, it suffices to consider Tr(£) and 
Tr\C)-Tr{C^). 

The following proposition tells us how to pose initial value problems for the 
Boussinesq system. The proof uses a different technique than the one used in 
[38], which is possible due to the fact that for the Boussinesq system initial 
values can be given on staircases. However, the staircases are not necessarily 
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standard staircases, which they would be in the framework of [38] (at least for 
equations defined on the square, such as the Boussinesq system). 

We call a staircase ascending, if it goes from the lower left to the upper 
right, that is, if it is a sequence of neighboring lattice sites with I and m non- 
decreasing. And we call a staircase descending if it goes from the upper left 
to the lower right, that is, if it is a sequence of neighboring lattice sites with I 
and —m nondecreasing. 

Proposition 3 The following initial values problems for the Boussinesq system 
are well-posed. 

• At every point on an ascending staircase take the components u,v of the 
vector u as initial values. 



• 



On a descending staircase take u, v, w at the lower left corners, v at the 
upper right corners, and u, v at the other points as initial values. 



Proof: The proof consist of two parts. Firstly, we show that the values at all 
points of the staircase can be obtained from the initial values. Secondly we show 
that any of the four vectors u, u, u, u, can be determined from the other three. 

i) For ascending staircases the first part is easy. Going along the staircase 
from the lower-left to the upper-right at each horizontal step the com- 
ponent w is obtained using equation (34o), whereas at the vertical steps 
equation (346) can be used. For a descending staircase we can do a similar 
thing, except at the upper-right corners. Equation (34d) can be solved for 
ui+i.m+i and used to get the u-components at the upper-right corners. 
Once u has been calculated w can be calculated in two ways, using either 
(34a) or (34b), leading to the same result. 

ii) The values of u can be obtained as follows. First calculate u from (34c). 
Then w can be obtained from (346) and v from the up-shifted consequence 
of (34a). We can obtain u in a similar way. This follows from the fact that 
interchanging the left-shift with the up-shift and 7 1— )• —7 is a (discrete) 
symmetry of the system. Finally, to obtain u one uses the consequence 
(34rf) to calculate u, after which w is found using a shifted version of either 
(34a) or (346), and v is calculated using (34c). Finally, due to the discrete 
symmetry which interchanges the up-shift with the down-shift and the 
left-shift with the right-shift, together with m -H- w and 7 h^ —7, it follows 
that u can be obtained from given values at the other sites. 

D 
It follows that, with s — (si, S2) G Z x Z such that S1S2 7^ 0, the dimension of 
an s-periodic reduction is 2(| si | + | S2 |). 
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(n-l,l)-reduction 



Xi X2 X3 Xi 

-m- ■ 



X2ri-2 
ff 




X5 xe 

. . .^. . . 



We take si — n — 1 positive (but 
n ^ 3), S2 = 1, and we consider the 
following initial value problem, with 
i = 0, 1, . . . , n — 1 and / G Z. 

U(n-i)i+id — a;2i+i, 

V{n-l)l+i,l ~ X2i+2, 



Figure 16: 2n initial values 



where the index on x is taken mod- 
ulo 2ri, see Figure 16. 
Using the equations (34a), (34c), and (34(i), in that order, we calculate, W(ri-i)(,i-i 

X2n-3X2n-l ~ X2n-2, M(„_l)i+l,i_l = Xi + 7P, and W(„_i);+i_;fc_i ^ X2 + X37P, 

where 

(35) 



1 



(36) 



a;2n-l(a;2n-3 - X3) - .T2„-2 + X^ 

The right-shift induces the 2n-dimensional mapping 02n :, 
Xi^Xi+2, i e {1,2, . . .,2n - 2}, 

X2n-1 1-^ Xi +"fP, 
X2n ^^ X2+ 7X3 P. 

The monodromy matrix is 

C = M(a;2„-i,w„,a;i,X2)L(a::2„-3,w„_i,a;2„-i,a::2„)- ■ ■ L{xi,wi,xz,Xi), 

where Wi = a;ia:;2„_i - a:;2„ and w^+i = a;2j-i.T2i+i - X2i, i = 1, 2, . . . , n - 1. 

The number of functionally indepen- 
dent integrals we have obtained for 
reductions with period s = (n— 1, 1) 
is given in table 2. From this table 
it seems we need to d-reduce map- 
ping (t)2n hy d — 2 dimensions, or, if 
3 divides n, by d = 6 dimensions. 



n 


2 


3 


4 


5 


6 


7 


8 


9 


# 


1 





3 


4 


3 


6 


7 


6 



Table 2: Number of functionally in- 
dependent integrals of (j)2n- 



The mapping (j)2n has two symmetries, generated by 

d ^ d 



Vl 



E 



dx' 



V2 



2i 



^-^ dx 



■a;2j-i 



2i-l 



dx 



2i 



The easiest way to check that these vector fields are generators of symmetries 
indeed is using the Jacobian, we have Jv = 4>2n{v) when v — (0, 1, 0, 1, ... , 0, 1) 
or f = (1, xi, 1, X3, . . . , 1, a;2n-i), where the Jacobian matrix of (/)2„ is given by 



J 



/2«-2 

h iP^H 
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where Ik is the k x k identity matrix and H is the 2 x (2n — 2) matrix 

X2n-1 -1 ••• -X2n-1 1 (x^ - X2n - 3) 

l/P + X3X2n-i -X3 ••• -a:;3a;2«-i x^ x^ix^ - X2n - i) 

The two symmetries of the mapping 02n correspond to the following symmetries 
of the original lattice system (34) : 

(u, V, w) H> (u, u + e, ui — e), (u, v, w) H> (?! + e, w + eu, li; + ew). (37) 

It can be verified that the functions 

Vi = X2i-i - X2i+i,yn+i-i = X2i - a;2i+2 +a;2j-i(a;2i+i - a;2j-i), 

with i — 1, . . . , n— 1, are joint invariants of the these symmetries and functionally 
independent. In the reduced variables y, we get a 2(n — l)-dimensional volume 
preserving mapping 

Vi^Vi+i, i e {l,2,...,2n-2},i^n-l, 

n-l 

yn-1 ^ -^yi-iQ, 

i=l 

n— 2 n—1 n—1 

y2n-2^ -C^yn+1 + yt+i X! yi)-iQ^y^-' 

i=0 j=i+l i=2 



where Q = J2z=i (yn+» + Vt+i J2j=t+i Vj) 



n — 3/ j^ Y^^^l 

i=l [Vn+i + J/i+1 Z^j=i+ 

When n — 2 the reduced mapping is, in terms oi X — — yi, 1" = — 2/2 

{X, Y) ^ {-X + j{Y - X^), -Y + X^) 

which carries the invariant y(X^ — Y) + jX, cf. [22, equation 5.31]. Note, in 
[22] the case si — $2 was studied, in particular the involutivity of the integrals 
was established in any dimension. The authors defer the actual counting of 
independent integrals to a future study. However, they also state that the 
investigation of lower-dimensional examples (with si = S2) indicate a sufficient 
number of invariants are functionally independent. 

We have verified that all integrals we found, see Table 2, survive the 2- 
reduction. Hence, in those cases, except when 3 divides n, the staircase method 
provides enough integrals for integrability. Next we will show that if 3 divides 
n, but n 7^ 3, we can further reduce the mapping by four dimension. 

The Boussinesq system (34) has some additional symmetries, which depend 
on the lattice variables (u — u; .„), 

(u, V, w) ^{u,v + eC'+"+\ w - eC'+"), 

(u, V, w)^{u + £('+", V + eC'+"+^w, w + eC'+""^u), (38) 
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where C is a primitive third root of unity, that is, C^ + C+ 1 = 0. The generators 
of the corresponding transformations acting on the initial values (16) are 





V3 


->;c 


dx2t' 


V4 = 




r 

V 


-1 5 

dx2i- 


- + CX2^- 
1 


d \ 

^ dx2i ) 


and. 


taking 


the conji 


ugate root C' 


> 












^5 = 


n 


d 

dx2i' 


V6 = 


71 


V 


+1 5 


- + C.^'X2^ 
-1 


d 




dx2i- 


'^ dx2i 



Now assume that 3 divides n. We construct real vector fields by taking the 
following linear combinations, with i = 1, 2, 3, 

n/3-l 

"/^^^ d d 

W3+t = iv2 + C^'W4 + Cva)/^ = E ^6j+2i-l 



.^Q dxQjj^2i dXej+2i+l 

where X2n+i = Xi. These vector fields are 3-symmetries of the mapping (j>2n- 
Let J^ be the Jacobian matrix of <j>2n- Note, the vector fields Wi, i ^ 1,4 can 
be obtained from Wi — J^(j)2n(wi-i)- According to [13, proposition 1] it suffices 
to verify that wi and W4 are 3-symmetries. Also note that J^(f)2nVi = Vi for 

ie{i,...,6}. 

The following polynomials form a complete set of joint invariants of the 
vector fields wi,W2, ■ ■ ■ ,we: 

Zi ^ X2i-i - X2i+5, I G {1, . . . ,n - 3}, 

Zn-3+i = X2i+l{x2i-l — 2:2^+5) — X2i + X2i+Q, i G {1, 2, 3}, 
Zn+i = {X2i-1 — X2i+5){x2i+6 — X2i+12) 

- {x2i -~ X2i+e){x2i+5 - X2i+ii), i G {1, 2, . . . , n - 6}. 

Another joint invariant is given by P, see (35), which therefore should be ex- 
pressible in terms of the Zi. However, we haven't found a general formula for 
P{z). In terms of the z- variables the mapping 0i2 reduces to 

Zi 1-^ Zi+i, i G {1,2,4,5} 

7 



Z3 1-^ -zi 



2^5 - 22^3 

Z2 



zq i~> -Zi + Z1Z2 -7 



Zh ~ Z2Z-i 
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and (l)2n, with 3 divides n > 4, 

Zi n- Zj+1, i G {1, 2, . . . , 2n — 7}, i j^ n — 3,i ^ n 

n/3-l 

i=l 

^K+1 + Z4{ziZ2 - Zn-2) 

Zn '^" J 

Zl 

Z2n-6 ^ A- -fBP, 

where A = {x2n-ii ~ X2n-b){x2n-A - X2) - {x2n-w " a:2n-4) (a;2n-5 " xi), and 
-B = a;3(x2n-ii — a;2n-5) — a;2„-io + X2n-i are also joint invariants of the six 
vector fields Wi. All the integrals we found survive the above reduction. Thus, 
as one can see in Table 2, also when 3 divides n the staircase method provides 
a sufficient number of integrals for the 6-reduced map to be integrable. 

We found two functions which are 3-integrals, that is, integrals of (/)2„. These 
are 

n/3 

h = 2j(^6:(+i - a;6j-5)(a:6j-2 - xqj-zxhj-i) + xqj-i{xqj+2 - 2:6^-4), 

n/3 
12 = 2J X(ij-i{xQj-^XQj-z - XGj-i - XGj+iXQj+3 + Xej+2) + XQj{xej+3 - Xqj-s). 

We have the following action of (j)2n on the 3-integrals 

(»i,«2) ^ (»2 - n +7,-n +7)1 (39) 

whose third power is the identity. Note that by applying the map to one of 
the 3-integrals gives us the other one but no third functionally independent 3- 
integral can be obtained in this way. The 3-integrals admit the vector fields Wi 
as symmetries and two of the n — 3 functionally independent integrals found by 
the staircase method can be written in terms of them. For n = 6 we have 

ii = -21^22:3 + Z1Z5 + Z3Z4, i2 = Z2Ze + Z3Z4. 

and two of the three functionally independent integrals found by the staircase 
method arc given by 

il + ii{ii -12- 7), *i(7 - «2)(7 - n + ^2) 

which are both integrals of the 2-dimensional map (39). An extra advantage of 
working with expansion (34) instead of traces of powers of C is that the third 
functionally independent integral factorizes nicely as 

2425^^6(^12:3 + ze){ziz^ - Z4Z2Z3 - Z2Z1Z5 + Z2Z1Z3 + az2)- 
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(n-l,-l)-reduction 



a;2n-2 




X4 X5 Xg 
V" 



Figure 17: Periodic initial value prob- 
lem for the Boussinesq system 



We take si = n — 1 positive (n 7^ 2), 
S2 = —1, and we consider the fol- 
lowing initial value problem, with i = 
0, l,...,n-2, fc e Z, 

U(n-l)k+i,-k = 2;2i+l, 

^(n-l)fc+i,-*; = 2;2i+2, 

^'(n-l)(fc+l),-fe = X2n-1, 

W(n-l)k-k = 2;2„, 



see Figure 17. 
We right-shift the initial values X2n-i and a;2n using equations (34(i) and (34a), 
respectively. Then we right-shift a;2„_3 a second time and are able to determine 
the image of a;2„_i, using equation (34c). Thus, we find the 2n-dimcnsional 
mapping C2«: 

Xi^Xi+2, i e {1,2,. . .,2n - 3}, 

a::2n-2 - X2 
a:2«-3 ^ , 

X2n-Z - Xi 

a:2n-2 ^ a;2n-i, 

(a;2n-3 - xi)(a;i(x2„-i - X4) - 7) 



a;2n-i 



a;2n 



-a;2n 



a;2n-2 - 2:2 - a;3(x2„-3 - Xi) 



-X2 + X1X3, 



which is measure preserving with density X2n-3 — xi (n ^ 2). The mapping C,2n 
admits the symmetries 



d 



dX2 



d 

dx2i 



E 



d 



dx2i 



and 



Xl- 1- C2n(a:2n-3)- 



9a;2„ 



9a;2n- 



n-1 

E( 



9a; : 



+ a;2i-i 



2i-l 



9a;2i 



Hence, it can be reduced to a (2n — 2)-dimensional mapping. We have verified 
up to n = 7, that the number of functionally independent integrals, obtained 
by fc-expansion of the coefficients in (34), with 

JO- = M^^{xi,X2n,Cnix2n^3),X2n-l)L{x2n-3,Wn~l,Cnix2n-3),X2n-l) 

n-1 

• Y\_ L{x2i-3,'Wi-.i,X2i-l,X2i), 
i=2 

is n— 1, except when 3 divides n+1, where the number is n — 3. Also we verified 
that all these integrals admit the above vector fields as their symmetries. If 3 
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divides n + 1 the symmetries (38) yield 3-symmetries of the mapping (^2n- As 
in the previous example we take linear combinations to get 



(n-2)/3 



d 



(n-5)/3-l 



dx, 



6j+2 



, W4^ { Y^ Xqj, 



d 



d 



j=o 



dx^ 



6i+2 



dx, 



X2n-3 



d 



6i+3 



dx2n-2 ' 



3=0 

together with Wi = J^C,2n(wi-i), i y^ 1,4. Taking n = 5 we obtain the following 
4-dimensional 6-reduced mapping, in terms of j/i — xj — xi,y2 = xq + xio — 
X5Xi,y3 ^ xs-X2+X3{-X7+xi),y4 ^ {-xr+xi){-XQ+X4-X5X3)+x<i{-xs+X2): 

(2/1,2/2,2/3,2/4) ^ —,2/3,— ,2/4 + 7 

V2/1 2/1 2/1 

This mapping is measure preserving with density y^. It has two 3-integrals, 
H = 2/22/3/2/1 s-nd 12 = 2/4, which satisfy 

(«1,«2) ^ («2,-«l -«2 -7), 

whose third power is the identity. The two functionally independent invariants 
found by the staircase method can be expressed in terms of the 3-integrals as 

i\ + i\+ 1x12 + 7(«i + «2), «i«2(n + «2 + 7)- 



4.3.1 The potential Korteweg-De Vries equation, (3,0)-reduction 

Consider initial values for the pKdV equation (22): uj.o — Xk with k = I mod 

3, k G {1,2,3}, as in Figure 18. 

Updating these by shifting them up- 
wards, and imposing the image to be 
periodic, we have to solve the system 

{Xi - 'X>2)(X2 -Si)= a. 
{x2''Xi)(xz-3>2) = a. (40) 

[xz -x{)(xy -xi) = a. 



Xi 


X2 


X3 


Xl 


X2 


X3 



Xl 



Xl 



Figure 18: Initial values on a 
horizontal staircase with periodicity 

for the image points xi, X2, x^. 
In terms of reduced variables qi ~ xi ^ X2, q2 ~ X2 — X3 the volume preserving 
two- valued correspondence is ip± : 



qi l-^' qi 

92 1-^ <?2 



a{a + q2{qi + 92)) + gig2(gi + g2)/2 + (gi/2 + q2)R 
{a + qiq2){a - qi{qi + ^2)) 

a{a + qi{qi + 92)) + gig|(gi + g2)/2 - (gi + g2/2)i? 
{a + qiq2){a - q2{qi + q2)) 



where 



R ^ ±v/4a3 + {qiq2{qi + q2)Y . 

This correspondence admits the integral qiq2{qi +92), which can be obtained 
by taking the trace of the monodromy matrix. 
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5 The staircase method on quad-graphs 

Recently, in [4], a geometric criterion was given for the well-posedness of ini- 
tial value problems on quad-graphs. In this section we show that for 'regular' 
quad-graphs, those that permit periodic solutions, the staircase method can be 
applied. We use equation H3s=o, from [2], 



J (a, 6, c, d) :~ p{ab + cd) — q{ac + bd) — 0, 



(41) 



on two different quad-graphs, cf. Figure 9d and 9e in [4]. 

We start with a brief introduction to the idea of a quad-graph. For a more 
thorough treatment and references to the literature we refer to [4] . A quad-graph 
is a planar graph with quadrilateral faces. Fields are assigned to the vertices 
and parameters to its edges. In the class considered in [4] opposite edges carry 
the same parameters and the (multi-linear) equation is supposed to have D4- 
symmetry (so that the equation can be defined on each face independently of 
its position in the quad-graph). Due to the first property there are sequences of 
adjacent quadrilaterals on which the value of the parameter is constant. These 
are called characteristics. The main result in [4] states that an initial value 
problem P is well-posed if and only if each characteristic intersects P in exactly 
one edge. 

An example of a well-posed initial value 
problem is given in Figure 19: from 
given values a, b, c, d one can calculate 
e, /, g successively, they are uniquely 
determined. Indeed, the three dotted 
lines are characteristics and they each 
intersect the initial value path in ex- 
actly one edge. 

In both Figures 20 and 21 a finite piece of a quad-graph is given. We as- 
sume these quad-graphs extend periodically in both the vertical and horizontal 
directions. For equations on such a doubly infinite quad-graph there exist a 
two-parameter family of periodic reductions, as for equations on the regular Z^- 
lattice. In fact, these quad-graphs are perturbations of the regular Z^-lattice, 
we call them irregular Zi^ -lattices. 

Theorem 1 extends to the more general setting of irregular Z^-lattices. The 
proof is similar as in the 1? setting: the fact the transfer matrix Ca.t does not 
depend on the path from a to b follows from the Lax-condition. 




Figure 19 
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In Figure 20 the lattice pa- 
rameters are attached to 
the edges as follows: p to 
the horizontal edges, q to 
the vertical edges, and r to 
the diagonal edges. From 
initial values a, 6, c, d the 
values e,f,g are deter- 
mined by 



Figure 20 



Qprid,a,e,b) 

Qqr{c,d,f,e) 

Qpq{e,b,f,g) 



0, 
0, 



We update values by shifting over (2, 1). Thus, the mapping is 

{a,b,c,d) 1-^ {c,f,g,b), 
where 



f = d- 



'•{rb — pd) + qc{rd — pb) 



rb{qa — pc) + qd{qc — ap) 



aq{+rb — pd) + rc{rd — pb) ' qb{qa — pc) + rd{qc — ap) 

A Lax pair L, M for equation (41) is given by L = La.b{p), M ~ La.dq), where 

1 f ka —pab 



La,b{p) = 



p —kb 



see [16]. The monodromy matrix Ld.a{p)Lc.d{q)Lb.c{p)La,b{'r) yields one inte- 
gral. 

We perform the following reduction: in variables x — a/ c, y ~ b/d the 
mapping is expressed 



{x,y) ^ 



p{qy + rx) — qiqxy + r) \ p{rx + qy) — r{q + rxy) 1 



^p{qx + ry) — q{rxy + q) x^ p{qx + ry) — r{qxy + r) y J 
and its integral is 

2 f^ , y 



P 



p{r + q) { x + y -\ \ — ] + { xy ^ ] qr. 

y x/ \ X yj \ xy^ 



In the next example we need four lattice parameters to ensure that the two 
lattice parameters on each quadrilateral differ. 
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Given initial values 
o, 6, c, d, the values 
e, /, g, h, i can be deter- 
mined by the following 
equations: 

Qrp{b,c,e,d) = 0, 

Qsp{a,b,f,e) ^ 0, 

Qrq{e,d,g,a) = 0, 

Qsq{f,e,h,g) = 0, 

Qrsig,a,h,i) = 0, 



QqpiaJ,j,h) 



0. 



The values a, b, c, d are repeated periodically on the 'staircase' which extends 
along a 'diagonal' of the quad-graph. The particular way of choosing lattice pa- 
rameters {q,p,p,q, see Figure 21, and we also take r,s,s,r vertically) ensures 
that the periodic solution has the same period as the initial values. It is im- 
portant to notice that when going one step to the right (on the 7? part of the 
quad-graph) the lattice parameters p, r are interchanged with q, s respectively. 
Therefore, we consider the mapping 

(a,b,c,d,p,q,r,s) h^ {h,i,a,j,q,p,s,r). 

After the transformation x — a/c, y — b/d we are left with 

[x,y,p,q,r,s) ^ ( , ,q,p,s,r), 

X Xd yyd 



with 



Xd 



Vd 



—qsp + pq{qs + pr)x + 2ypsrq — [rq + ps) (qs + pr)yx — y sr q 

22, / I \2|/22|2 2\2 222 

—pq X r + sryqs + pr)y x + (p q + s r )yx — y ps x r, 

y rqp x — y rq p — y qsx p + y q sxp + yr s — ypqr x + 2yrqsx p 

—yp sxr — yq sxr + yq p — yps xq — r qsx + r spx + rs xq — ps r, 

(s p — s{qs + pr)x — s ry + (prq + s )yx + x rqs — yqsx p) 

■{—qsp + q{qs -\- pr)x + ysrq — (srp + q )yx — q x r + x q py), 

{—qsp + {prq + s^)x + ysrq — s(qs + pr)yx — x'^s'^r + x'^s'^py) 

■{q p — {srp + q )x ~ yq r + q{qs + pr)yx + x rqs — yqsx p). 



The monodromy matrix, 

Ld.a{l)Lc.d{p)Lb.c{^)La,b{s) 

yields one invariant, in reduced variables, 

x{l + y'^){qs + pr) + y(x^ + \){pq + sr) — {x^ + y'^)rq — {x'^y'^ + l)ps. 
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The linear part of the mapping is easily solved by 

P - Ci + C2(-l)", q = c,- C2(-l)", r = C3 + C4(-l)", s = Cg - C4(-l)" 

where n is an integer counting the iterations of the map. Thus we may obtain 
an alternating two-dimensional map with one integral, cf. [28]. However, this 
map would take too much space to write explicitly. For special values of the 
parameters and one of the initial values we have observed that the growth of 
this map is quadratic, which indicates it is integrable. 

6 Multivaluedness of iterates of correspondences 

We have seen two examples in which correspondences arose: (n, 0)-reductions 
of the QD-system in section 3.2.1 and of the pKdV-equation in section 4.3.1. 
In general the multi-valuedness of the iterates of a correspondence would grow 
exponentially. However, for integrable correspondences one expects the multi- 
valuedness to grow polynomially instead. 

6.1 The pKdV (3,0)-correspondence 

Generically the multi-valuedness of the 
/th iterate of a two-valued map would 
be 2K This is not the case here. Due to 
the X o 2/ symmetry of the system (40) 
the correspondence (p equals its own in- 
verse. The relations ip^ = ip^ imply 
that the /th iterate of the correspon- 

Figure 22: Polynomial growth of deuce is (/ + l)-valued, see Figure 22. 

multivaluedness of iterates. All points on the same vertical line have 

the same value. 

6.2 The QD (n, 0)-correspondence 

As one can easily verify the mappings r, cr, given by equations (20,21), satisfy 

TdT ~ aacr and raa = aar. (42) 

Due to these relations the Ith iterate of the correspondence is 21 valued, see 
figure 23. 

Proposition 4 With compositions of mappings (42), the Ith iterate oj the cor- 
respondence (cr, t) is 21 valued. 

Proof: The relations (42) can be used to rewrite every word in two symbols a 
and T of length / into either a word that docs not contain ra or a word that 
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TTTTO 




TTTTT 



Figure 23: Graphical representation of Proposition 4 



only contains tu at the end. There are 21 such words: 



TTT ■ ■ ■ TTT, aTT ■ ■ ■ TTT, GGT ■ ■ ■ TTT, . 



, aaa ■ ■ ■ aar, 



TTT ■ 



■ TTG, (TTT ■ 



■ TTa, auT ■ ■ ■ TTa, . 



, aaa ■ ■ ■ aaa. 



So there are at most 21 incquivalent words of length I. To show that there are 
exactly 21 inequivalent words of length I we proceed by induction. We suppose 
there are are exactly 21 inequivalent words of length I. For any three words 
u, V, w we have uv = uw ^ v ~ w. So if two different words of length I extend 
to an equivalent word of length Z + 1, this word can be written (we concatenate 
from the left) as, u — tv = aw. If a third word of length I would extend to u 
then we have u — tz or u ~ az which would imply z — v or z — w respectively. 
Therefore words of length I that extend to words of length I + 1 coincide at most 
pairwise, giving an lower bound of 21 on the number of words of length I + I. 
However, since there are two words of length / + 1 that are not equivalent to 
any other word, namely ttt ■ ■ ■ ttt and aTT ■ ■ ■ ttt, there are at least 21 + 2 
inequivalent words of length I + 1. D 



The same relations hold in the 2n-dimcnsional case. 

Proposition 5 The 2n- dimensional mappings t — T2n (17) and a = a2n (13) 
satisfy the relations (i2). 

Proof: That the second of the relations (42) holds is easily established as ctct, 
i.e. the mapping Xk ^^ Xk+2, clearly commutes with r. 

For the first relation we solve the equations (12) for the Xi in terms of the 
Xk to find the inverse t~^, with i = 1, 2, . . . , n. 



a;2j-i H^ a;2j-i 



W^i 






X2i H> X2i 



^W^^-2' 
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where 

^fe" ^ E n ^2.-i+/c n ^2.+2+fe . (43) 

One can verify that <j{W^) = Q^. This implies that ctt~^ is given by 

V2i-2 

a;2i i-^ a;2i-3 7^;i — , i = l,...,n 

which equals its inverse tct^^. Multiplying rcr^^ = ctt^^ from the left by ra 
and from the right by cr, using rcrcr = crcrr, gives us rar — aaa. D 

Remark: One can also directly prove that the t^^ provided is the inverse 
of T. This relies on the fact that t{W^) = QJJ, which in terms of Q amounts to 
the identity 

n — l n i—1 n—1 n—2 

n ^2. = E(n ^2,-i)(n ^2,+4)(n q2,+2.), 

which can be proved using 

n— m+l 
i=3 

which in turn relies on a generalization of the relations (16), 

m+l n—l 

Qr" - Q^„"( n ^2.) + (Q^ n ^2.+i). 



7 Concluding remarks 

We have shown that the staircase method provides integrals for mappings and 
correspondences derived as s-periodic reductions of lattice equations and sys- 
tems of lattice equations. We also showed that such mappings and correspon- 
dences can be order- reduced systematically, using symmetries of the lattice equa- 
tions. In all examples encountered the staircase method yields sufficiently many 
functionally independent integrals for the d-reduced mappings and correspon- 
dences to be completely integrable. However, we know the above statement is 
not true in general. In [34] periodic reductions of systems of PAEs are consid- 
ered for which the staircase method does not provide sufficiently many integrals. 
However, in those cases it was observed that a certain linear combination of in- 
tegrals factorises into a product of 2-integrals, from which then another integral 
can be obtained. 
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For posing initial value problems we used the method laid out in [38]. Al- 
though, for the Boussinesq system we presented an alternative approach. We 
have calculated integrals for mappings/correspondences with dimension up to 
say 20. Therefore we were able to check functional independence of the integrals 
using the symbolic software package Maple [18]. 

Closed form expressions for integrals of mappings with arbitrary dimensions 
(corresponding to reductions with S2 = —1) have been given in [40, 36]. There, 
the lattice equations are of Adler-Bobcnko-Suris type [2, 3] and the integrals 
are expressed in terms of multisums of products. Their functional independence 
and involutivity is being investigated [37]. 

It is an open problem whether such closed form expressions can be obtained 
in general. In particular, we don't know whether such expressions can be given 
for the Bruschi-Calogero-Droghei equation, the Quotient-Difference algorithm 
and the Boussinesq system presented here, except for QD in the case si = 
0. Another question is how to obtain symplectic structures for the mappings 
studied in this paper. This would enable one to conclude complete intcgrability. 

We have obtained a few results on mappings of arbitrary dimension. These 
include an explicit (2 + (— l)")-reduction of the mapping related to the (n— 1, 1)- 
reduction of the lattice potential KdV equation, as well as an explicit expression 
for a 2-integral of the mapping in the case that n is even. We proved that the 
(n — 1, l)-reduction of the Boussinesq system can be d-reduced with d = 6 if 
3 divides n and with d — 2 otherwise. Also, we presented two 3-integrals for 
this mapping (expressed in terms of the original variables). For the (n — 1, — 1)- 
reduction of the Boussinesq system we have showed that the mapping can be d- 
reduced with d = 6 if 3 divides n+1 and with d = 2 otherwise. Finally, the (n, 0)- 
reduction of the QD-system yields a 2n-dimensional 2-valued correspondence. 
We have given an explicit expression for this correspondence and showed that 
its nth iterate is 2n-valued. It would be interesting to investigate other ways of 
establishing such a result because one does not always have explicit expressions 
at hand. 
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